function plotdom3(A)

M = max(A);
A(:, 3) = -A(:, 3);
M(3) = -M(3);

A = sortrows(A, -3);

colormap(gray)
caxis([0 1])

domdown([], A, M);
view(-45, 30);
set(gca, 'XGrid', 'on');
set(gca, 'YGrid', 'on');
set(gca, 'ZGrid', 'on');

hold on
scatter3(A(:, 1), A(:, 2), A(:, 3), 40, 'r', 'filled');
hold off

title('Pareto-dominated volume')
xlabel('Processing time (s)');
ylabel('Area (mm^2)');
zlabel('-Energy (J)');
